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j^ ' 1. Introduction. In this note we will study a special class of solutions of the three- 

O ■ dimensional steady-state Navier-Stokes equations 

(N 

S^: NSE 

<. divw = 0. ^ ^ 

! The equations have a non-trivial scaling symmetry u{x) — > Xu{Xx) and it is natural to try 

jr~i| to find solutions which are invariant under this scaling. Explicit examples of such solutions 

^ I were first calculated by L.D. Landau in 1944 ([L]) and can be found in standard textbooks. 

(See, for example, [LL], p. 82, or [B], p. 206.) The main idea of Landau's calculation 

is that if we impose an additional symmetry requirement, namely that the solutions are 

axi-symmetric, the system (NSE) reduces to a system of ODEs which, surprisingly, can be 

solved explicitly in terms of elementary functions. The solutions were also independently 

found by H.B. Squire in 1951 ([Sq]). More recently, the topic has been re- visited in [TX] and 

'i> ', [CK], where issues concerning Landau's solutions are addressed from a slightly different 

\^ I viewpoint. 

vn ■ 

T^ ■ In this note we prove that even if we drop the requirement of axi-symmetry. Landau's 

solutions are still the only solutions of (NSE) which are invariant under the natural scaling. 
Q ■ More precisely, we will prove the following: 

'^ Theorem 1. Assume that u: R^ \ {0} -^ R^ is a non-trivial smooth solution of (NSE) 

satisfying Xu{Xx) = u{x) for each A > 0. Then u is a Landau solution. In other words, 
u is axi-symmetric and, in a suitable coordinate frame, is described by formulae (E7) in 
Section 4. 

X' 

;-H I The proof of the theorem shows a connection between the scale-invariant solutions of (NSE) 

and the conformal geometry of the two-dimensional sphere. In fact, once the connection is 
understood, the formulae for Landau's solutions can be derived without much calculation, 
using just the geometrical properties of the two-dimensional sphere. 
Some implications of Theorem 1 are considered in the next two sections. 

2. Implications for regularity of very weak solutions. By a very weak solution of the 
steady-state Navier-Stokes system (NSE) in a domain Q G R^ we mean a divergence- free 
vector field u = {ui, U2, U3) G L^^^(O) which satisfies 
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/ [Ui^ipi -f UiU 

Jo. 



'' dxi 



1 



for each smooth, compactly supported, divergence-free vector field </? = ((/?i, Lp2i 9^3) in ^■ 

It seems to be an open problem whether very weak solutions of (NSE) are regular. Standard 
regularity theory can be used to show that very weak solutions are regular under the 
additional requirement that u G L^^^. The equations (NSE) are usually considered with 
the assumption that Vw G L^^^, in which case regularity follows by a straightforward 
bootstrapping argument. The assumption Vti G L^^^ is of course very natural when 
considering solutions describing real physical flows. However, one can speculate that very 
weak solutions might arise from a blow-up procedure of the usual weak solutions of the 
time-dependent Navier-Stokes equations at a possible singularity (if a singularity exists). 
The time-dependent 3-dimensional Navier-Stokes equations are supercritical with respect 
to the natural energy estimates, and in a blow-up procedure the information about energy 
can be lost. 

A natural first step in understanding the regularity of the very weak solutions above is to 
study the scale-invariant solutions in B? which are smooth in i?^ \ 0. A calculation (which 
can be found in [B], p. 209, and also in [T], and [CK]) shows that Landau's solutions 
are not very weak solutions of (NSE) across the origin. Therefore Theorem 1 implies the 
following. 

Corollary. Let u he a {—1) -homogeneous very weak solution of the Navier-Stokes equa- 
tions in R^, which is smooth away from the origin. Then u = 0. 

This result rules out only the simplest conceivable singularity of a very weak solution. For 
example, the question if one can have a non-trivial very weak solution smooth away from 
the origin and satisfying |w(x)| < C|x|~^ in R^ is not answered by Theorem 1 and - as far 
as I know - remains open. 

Remark: The question about the existence of (—1) -homogeneous solutions of the system 
(NSE) in R^ \ {0} (smooth away from the origin) can also be posed for n > 4. It turns 
out that for n > 4 there are no such solutions. This was proved independently by several 
authors ([St], [Sv], [T]) and it also follows from results in [FR]. 

3. Asymptotic behavior of solutions in exterior domains. Theorem 1 has some 
relevance for the problem of long-range behavior of solutions of the Navier-Stokes equations 
in exterior domains. (See for example [G] for an overview of this topic.) For simplicity, we 
consider here only the following special case. Let / be a compactly supported vector field 
in R^ and consider the equations 

-Au + uVu + Vp = f , 

in R-^ , 
divu = 

together with a "boundary condition" at 00, which might take the form u{x) -^ at 00 
and /^3 |Vwp < 00, for example. The existence of such solution was proved in a classic 
paper by Leray ([Le]), but there are many open questions about the behavior of these 
solutions for large x, see [G]. Theorem 1 implies, roughly speaking, the following: 



If a solution of the above exterior problem is asymptotically {~1) -homogeneous, then the 
terms of order \x\~^ must be given by a Landau solution. 

To give this a more precise meaning, let us consider the scaled functions ux and f\ defined 
by ux{x) = Xu{Xx) and fx{x) = X^f{Xx). The functions ux and fx satisfy the same 
equations as u and /. Moreover the functions fx converge to a distribution /, given 
by f{x) = M(x), where b = fj^^ f and d is the Dirac function. Assume now that ux 
converges to a limit u in, say, Lf^^{R^ \ {0}) as A ^ oo. Our notion of "asymptotically 
(—1)- homogeneous" used above can be definedhy requiring that this is really the case. (Of 
course, it is a difficult open problem to decide whether, for general large data, this might 
always be the case.) The limit functions u and / will again satisfy the same equations (in 
the sense of distributions). Under our assumptions the function u is smooth away from 
the origin, satisfies Xu{Xx) = u{x) for each A > 0, and, by Theorem 1, must therefore be 
a Landau solution or vanish identically. (The direction of the vector b will be the axis of 
symmetry of the solution.) For 6 = we will have -u = 0, which means that, under the 
above assumptions, the solution u decays faster than \x\~^. 

4. Proof of Theorem 1. Let m be a (—1) -homogeneous vector field in i?"^, smooth away 
from the origin. Clearly u is determined by its restriction to the unit sphere S'^ C R'^. 
For X E S'^ we decompose u{x) as u{x) = v{x) + f{x)n{x), where n{x) = x is the outer 
unit normal to 5"^, and v{x) is tangent to S'^ at x, i. e. v{x) ■ n{x) = 0. We now write 
down the Navier-Stokes equations for u and as a system of PDEs on S"^. If tt satisfies the 
Navier-Stokes equation in R^ \ {0} in the very weak sense defined above, it is easy to see 
that there exists a suitable pressure function p in R'^ \ {0} which is (— 2)-homogeneous 
and smooth away from the origin. The function p is also determined by its values on S'^, 
and the system (NSE) can be written down as a system of PDEs on 5"^ for v,f and p. 
The differential operators in what follows will all be differential operators on 5"^, defined 
by the usual conventions of Riemannian geometry. The differential forms on S"^ will be 
identified with vector fields and vice-versa, as is usual on Riemannian manifolds. The 
Hodge Laplacian dd* + d*d on 1-forms will be denoted by —Ah. (The reason for writing 
it as —Ah, with the minus sign, is to keep the equations on 5"^ in a form which resembles 
the standard euclidean form of the equations as much as possible.) The Navier-Stokes 
equations (NSE) for u written in terms of v, f and p as equations on 5"^ are as follows: 

-A^v + vVv + V(p - 2/) = , 
-A/ + W/-/2-|t;|2-2p = 0, (S) 

divt' + Z = 0. 

A straightforward (although perhaps not the most illuminating) way to derive these equa- 
tions is to write the system (NSE) in spherical coordinates (see, for example, [B], p. 601) 
and check that for (— l)-homogeneous vector fields it reduces to the system (S). We remark 
that the spherical coordinates version of (NSE) in the second edition of the book [LL] (p. 
49) contains a misprint in the right-hand side of the first equation, where an incorrect 
expression sin^ 9 appears instead of the correct sin 6*. 



We will denote by uj the function on 5"^ given by dv = ujQq, where Oq is the canonical 
volume form of 5"^. This corresponds to the formula u = curlz; used in R^. 

By taking d of the first equation of the system (S) we obtain 

~Auj + div (vuj) =0. (El) 

Lemma 1. With the notation introduced above, we have a; = 0. 

Proof. Let L be the differential operator defined by Lw = —Aw + dw{vw). The adjoint 
operator L* is given by L*w = —Aw—vVw. The kernel of L* consists of constant functions, 
as can be seen from the strong maximum principle. The kernel of L must therefore also 
be one dimensional. Let us denote by wq a non-trivial function in the kernel of L. If 
wq changed sign on 5^, we could find a strictly positive smooth function h on S'^ with 
fg2 woh — 0. But this would mean that the equation L*wi = h has a solution. However, 
the last equation cannot be satisfied at points where wi attains its minimum. From this 
we see that the function uj cannot change sign. At the same time, the definition of lv 
immediately implies that J^a cu = 0, and we see that to must vanish. 

Remark. I assume the above argument is known in one form or another, but I was not 
able to find a good reference for it. 

Once we know that dv = 0, the first equation of (S) simplifies. Indeed, when dv = we 
have —A^v = — Vdivz; = V/, and we also have vVv = V|i'p/2. Using this, the first 
equation of (S) implies 

-|t'P+p-/ = c, 

where c is a constant. The second equation of (S) now gives 

-A/-2/ + div(/t;)=2c. (E2) 

Integrating (E2) over 5"^ and using the third equation of (S) we see that c = 0. Since 
dv = we can write v = Vf for a suitable smooth function (f on S'^. The equation (E2), 
together with the third equation of (S) and the fact that c = now gives 

A2(/; + 2Aifi - div ( AvjVv?) = . (E3) 

Letting w = 2 — A(f, the last equation can be re-written as 

—Aw + div (S/i^w) = . 

The solutions of this equation are well-known: They are functions of the form cie'^, where 
Ci is a constant. (An easy way to verify this is for example the following: Write w in 
the form Ci{x)e'-'^^^K We get an equation for ci for which the strong maximum principle 
implies that the solutions are exactly ci{x) = const.) Integrating w over the sphere we see 
that ci > 0. Hence we have 

-Av7 + 2 = cie'^ 



for a constant Ci > 0. Changing </? by a constant, if necessary, we can assume ci = 2 
without loss of generahty, and we end up with 

-Av? + 2 = 2e'^ . (E4) 

The interpretation of equation (E4) is weU- known (see, for example, [CY]): Let g be the 
canonical metric on 5"^ and let g be the metric on S"^ defined hj g = e^g. Equation 
(E4) says exactly that the Gauss curvature of the metric (7 is 1, i. e. the metric g is 
isometric to the metric g. In other words, we have g = h*g (puUback of g by h) for 
a suitable diffeomorphism h oi S"^ . From the definitions we also see that h has to be 
conformal or anti-conformal. Anti-conformal maps can be obtained from conformal maps 
by a composition with an isometry, and hence we can only consider the case when h is 
conformal. For a given conformal /i, the function ip is given by 

V7(a:) = log|/^'(x)|^ (E5) 

where h'{x) denotes the (complex) derivative of h at x. It is well-known (see e. g. [DFN]) 
that all conformal diffeomorphisms of 5"^ can be produced as follows. Let P:S^ ^ C 
be the standard stereographic projection, and let My. C ^ C be defined by 2 ^ \z. 
Let h\ = P~^ o Mx o P. Then any conformal diffeomorphism of 5"^ can be produced by 
composing a suitable h\ (with A > 0) with isometries of S'^. If </? is given by (E5) and we 
compose h with and isometry, then the function </? either does not change or only changes 
by being shifted by the isometry. Therefore in a suitable coordinate frame all solutions <f 
of (E4) look like the solutions generated by the special hx above. We now consider the 
standard spherical coordinates (6*, -i/;) on S"^, given by 

xi = sin 6* cos V' , 

X2 = sin 6* sin -!/; , (SC) 

xs = cos 9 . 

We will use the usual notation 60 = ^ for the tangent vector field on S'^ corresponding to 
^. Letting A = e~'^, calculating the maps hx above in these coordinates, and using the 
formula (E5), we obtain 



'^{x) = —2 log (cosh K — sinh k cos 9) . (E6) 



This gives 



d(p —2 sin 6* 

89 coth K — cos 9 

f = -A^ = 2e'^ -2 



(E7) 



(cosh K — sinh k cos 6*)^ 

which agrees with the formulae in [B], p. 207 if we set coth k = 1 + c and with the formulae 
in [LL], p. 82, if we set coth^ = A. The proof of Theorem 1 is finished. 



Remarks: 

1. As we already mentioned in Section 1, the Landau solutions (given by (E7)) do not 
satisfy the Navier-Stokes equations (NSE) across the origin. A calculation in [B], p. 209, 
shows that for Landau's solutions we have, in distributions, 

—Aw + div (u^u) + Vp = bS , 

where d is the Dirac function and h = h^n) is a non-zero vector in B? depending in a non- 
trivial way on the parameter n which parametrizes the solutions in the above coordinate 
frame. The exact formula for h can be found in [B], p. 209, and was also calculated in 

[CK]. 

2. li h: S"^ ^ S"^ is a non-trivial holomorphic map (which is not necessarily a diffeo- 
morphisms) the formula (E5) gives a function ip which is regular away from a finite set 
tti, . . . , ttm £ S"^ where h' vanishes. The function ip will generate a (—1) -homogeneous so- 
lution of the Navier-Stokes equations in the region R^ \ (U^Z^i?-(- ■ Oj), where Rj^ = [0, cxd). 
However, the vector field will not be locally square integrable in R^, except for the case of 
Landau's solutions, when h' does not vanish at any point. 



5. Open problems. An interesting problem is to try to repeat the above analysis when 
R^ is replaced by the half-space R\ = {x G -R"^, X3 > 0} and the boundary condition 
-u = is imposed on dR\ \ {0}. A simple calculation shows that there are no non-trivial 
axi-symmetric ( — 1) -homogeneous solutions in that case. However, it seems to be an open 
problem if this conclusion is still true without assuming the rotational symmetry. We 
refer the reader to the very interesting paper [Se], where a related situation is studied in 
a different context. 

Another interesting question is the following: Among smooth vector fields in R^ \ {0} 
satisfying \u{x)\ < C\x\~^ for some C > 0, are the Landau solutions the only ones which 
satisfy the Navier-Stokes equations (NSE) in i?^ \ {0}? Such questions are relevant for the 
problem of asymptotic behavior of steady-state solutions in exterior domains mentioned in 
Section 3. A first natural step in addressing this question is to look at possible infinitesimal 
deformations of Landau solutions in the above class. This leads to linear equations which 
can be reduced to ODEs by classical methods of separation of variables, due to the symme- 
tries of Landau's solutions. Based on numerical experiments with these ODEs, the author 
conjectures that the Landau solutions are rigid with respect to infinitesimal deformations, 
i. e. it seems that there are no new solutions bifurcating from Landau's solutions. 
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